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Universal scaling behavior of the single electron box in the strong tunneling limit 
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We perform a numerical analysis of recently proposed scaling functions for the single electron 
box. Specifically, we study the "magnetic" susceptibility as a function of tunneling conductance 
and gate charge, and the effective charging energy at zero gate charge as a function of tunnehng 
conductance in the strong tunnehng hmit. Our Monte Carlo results confirm the accuracy of the 
theoretical predictions. 



I. INTRODUCTION 

A single electron box is a low-capacitance metallic is- 
land, connected to an outside lead by a tunnel junction. 
It has first been experimentally realized by Lafarge and 
co-workers [H . Such a device exhibits Coulomb-blockade 
phenomena due to the large charging energy, which in- 
fluences single electron tunneling. By applying a gate 
voltage Vg it is possible to induce a continuous polar- 
ization charge na = CaVa/e on the island (see Fig.[T]). 
For a very small dimensionless tunneling conductance 
a = h / {2tt'^ Rt) ^ 1 [Rt is the resistance of the con- 
tact) the energy spectrum is given by the set of parabo- 
las Eninc) = Ec{n — uqY , where the integer n denotes 
the number of excess electrons. Ec = e^/2(Ct -f Gg) is 
the single electron charging energy, which sets the energy 
scale. As the tunneling conductance is increased, quan- 
tum fluctuations lead to a renormalized effective charging 
energy E'^. The properties of the electron box for a < 0.8 
are well described by perturbation theory in a, and these 
results have been verifled numerically ■ 
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FIG. 1: Circuit diagram of the single electron box. The box 
with excess charge n is indicated by the dashed rectangle. It 
is connected to a voltage source through a capacitor Co and 
a tunnel junction with resistance Rt and capacitance Ct- 

Over the last decade, the analytical expressions de- 
scribing a single electron box in the limit of large tunnel- 
ing conductance have been the subject of controversial 



debate Numerical investigations which were 

performed to test the numerous incompatible predictions 
have proven difficult or inconclusive, and in some cases 
even added to the confusion [sllsllgl lTollllI . It is the pur- 
pose of this paper to remedy this unsatisfactory state of 
affairs and present a comprehensive study of recently ob- 
tained scaling formulas, using state-of-the-art numerical 
techniques. 



II. THEORY 

We address the limit of large tunneling conductance in 
the framework of the Ambegaokar-Eckern-Schon model 
[T^fisj . The partition function for the model can be writ- 
ten as a path integral over an angular variable (f) (conju- 
gate to the number of excess charges on the island) 



Z = j V4) e-^I-^l , (1) 

where the Matsubara effective action reads {h= \ and j3 
denotes the inverse temperature) 
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The first two terms account for the charging effect, 
whereas the third, "dissipative" term describes the elec- 
tron tunneling. Introducing a two component unit vector 
S = (cos (j) , sin (j)) at any point of the Matsubara circle, 
A[(l)] can be viewed as the classical energy of the XY 
model with long range interaction 14]. The first term in 
Eq. (m is then identical to the spin wave approximation 
to an extra short range interaction. Finally, the term 
cx UG acts like a purely imaginary external torque on the 
XY spins. 
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We are primarily interested in the "magnetic" suscep- 
tibility and effective charging energy of the model 
By magnetic susceptibility we mean the quantity 



Xr, 



^0 



drdr' (5,(t)5,(t')) 



(3) 



characterizing the linear response of the spin chain to an 
external "magnetic" field coupled to Sx = con cj). The 
effective charging energy, in turn, is defined as 



(4) 



In the terminology of QFT a quantity like ^ is usually 
referred to as a topological susceptibility because of the 
topological nature of the second term in Eq. ^ , and the 
parameter 



9 = 2-imc 



(5) 



is called the topological angle. 

When both f3Ec , a ^ oo the model ^ develops a uni- 
versal scaling behavior characterized by the energy scale 
E* (Kondo temperature). In this regime = (27r^a)~^ 
plays the role of the bare perturbative coupling and the 
first term in Eq. ([2]) just provides an explicit ultraviolet 
cut-off of the dissipative action, with the cut-off energy 



A = 27r2 Ec a 



(6) 



The Kondo temperature admits a perturbative expansion 
of the form [Hi] 



E* ~ 2tt^Ec e-°"'' U + Y.^'^ 



k=l 



(7) 



where the symbol ~ stands for asymptotic equality. 
While the overall factor 27r^ in Eq. ([7]) is a matter of con- 
vention, all the expansion coefficients bk are determined 
unambiguously within the standard perturbation theory. 
In particular 
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(8) 



In the scaling limit, i.e., in the limit when A — > oo, 
a — *■ oo with the energy scale E* kept fixed, the mag- 
netic susceptibility Q is expected to possess the normal 
scaling behavior of the form Q 



1 



Xr, 



27r2a 



i^„(K,0) + O((/3A)-i) 



(9) 



where the parameter k is the inverse temperature mea- 
sured in units of the Kondo temperature, 

K = f3E* , (10) 

and the topological angle 9 is given by Eq. 



The scaling behavior of the effective charging en- 
ergy dH) is a more delicate issue. Let us note at this point 
that field theories possessing topologically nontrivial clas- 
sical Euclidean solutions (instantons) usually suffer from 
specific "instanton divergences" . The phenomenon was 
orig inally observed in 4D QCD [H]. Later it was exten- 
sively studied in the context of the 2D 0(3) nonlinear 
cr-model [16,17,18,19]. For Ec = oo, the model (HD ad- 
mits instanton solutions which result in the "small 
instanton problem" (see the appendix below). For large 
but finite Ec the Coulomb term not only makes the the- 
ory perturbatively well-defined, but also regularizes the 
"small instanton divergence". As Ec oo, the diver- 
gence restores and entails the anomalous scaling behavior 
for the effective charging energy 



E^, 2Tr^E* 



L cos{9) + Ft{K,9) +o{l)' 



(11) 



Here L is a temperature- independent constant, diverging 
when a — *■ cxD and o(l) is a correction which vanishes 
in this limit. The scaling function Ft is defined mod- 
ulo a transformation Ft ^ Ft + c cos 9, since any finite 
constant c can be absorbed by the divergent term. This 
ambiguity can be resolved by imposing a proper normal- 
ization condition, say. 



lim Ft{K,0) = 



(12) 



As soon as the normalization of Ft is chosen, the di- 
vergent term in Eq. (jlip is also defined unambiguously. 
According to the result of Ref. [7] 

L{a) = 2n'^a- 5 loga- C . (13) 

The constant C has been found recently in Ref. [2^ : 

C = 57B + 6 log2 + 10 logTT 18.492.. . , (14) 

where 7^ = 0.5772 ... is Euler's constant. 

As was argued in Refs. 7l l2lll23 |. the scaling behavior 
of thermodynamic functions of the model Q can be de- 
scribed in terms of solutions of the differential equation 

- 92 + ^2 exp { ey)-K^ sm^{9) 1 ^(y) = . (15) 

Namely, let 5'+ and be solutions of Eq. ([T5)l which 
are fixed by the asymptotic conditions 



*-(y) 



and 



TT 

— exp 

K 



^yK cos 9 



ie^ 



y 



«Ei(V) 



-00 



(16) 



(17) 



with Ei(z) — Then the scaling function Ft 

defined in Eq. (|lip , satisfying the normalization condition 
(HH), reads as follows 



Ftiii, 



(18) 
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where 

and W^[*+, = '^+dy'^- - is the Wronskian 

of the solutions '^+{y) and \E'_(y). It should be empha- 
sized that ^_ can be defined by means of the asymptotic 
condition (dH) for < < f only. For | < 6* < tt, the 
solution grows at large negative y and the asymp- 
totic formula ([TB)) does not define VP- unambiguously. It 
is possible to show that the function \E'_/r(l -I- 2k cos 61) 
is an entire function of the complex variable ^ = cos 9 for 
real y. So the solution '^-{y) can be introduced within 
J < 9 < TT through analytic continuation with respect to 
the variable 9 from the domain < 9 < ^. 

The scaling functions Fjn(K, 9) defined in Eq. ([3 were 
proposed for the cases = and 6* = tt in Refs. [3| and 
[2l|, respectively. These results can be summarized as 
follows. Let us introduce the function 



M{K,9)=A+{y)+A^{y) + 



dyA+{y) - dyA^jy) 
s+{y) - s-{y) 



(20) 



where s±{y) = 9ylog^±(y) and 



du 



-oo J- 



dwe"*i(w) 



(21) 



Notice that the RHS in Eq. ([20|) does not actually depend 
on the choice of y. Then 

F,niK,9) ^ M{k,9) with 61 = and n . (22) 

It is essential to point out that no arguments have been 
given that relation ([221) holds true for any values of 9 
different from 9 = and 9 = tt. 

Let us now turn to the expression for the effective 
charging energy at zero temperature and gate charge. 
Eqs. (|lim3p imply that 



Ec 
where 



/(«) 



oo 



(23) 



f{a) = A^^a^ [ 1 - 1, + 0[a-^)] [ L{a) + o(l) ] , (24) 



8a 

and L{a) is given by Eq. (fTS]) . 

There is a controversy in the literature about the 
asymptotic behavior of E}j at no = = in the large- 
a limit. All the predictions agree on the exponential sup- 
pression of the effective charging energy, but there are 
numerous conflicting results for the leading term of the 
pre-exponential factor /(a), with powers ranging from 
a to a Mil. Several attempts to resolve the is- 
sue by numerical simulations [sUslllOjl have failed because 
of the difficulty to reach the large-a regime where the 



exponential factor in Eq. ((23)1 dominates. Using efficient 
algorithms, it is possible to compute the effective charg- 
ing energy at essentially zero temperature for a < 2 [ll| . 
In this range of a's, there was no obvious agreement with 
any of the theoretically predicted asymptotic formulas. 
While the pre-exponential factors of Refs. dUIil could be 
ruled out on the basis of the numerical results presented 
in Ref. [l3| , no definite conclusion could be reached con- 
cerning the predictions of Refs. [1|7|- It should be noted 
that the value of the constant C defined in Eq. (fH)) was 
not available at that time. 

The goal of the present work is to verify the analytical 
results quoted above, including expression (I24|) . by means 
of Monte Carlo simulations. 



III. NUMERICAL METHOD 

In order to simulate the electron box, we map the sys- 
tem ^ for no = to a finite chain of classical XY 
spins, by discretizing imaginary time into N equal slices 
Ar = /3/7V, 



1 



N 



2Ec^r 



0fe)) 



an 



E 

k<k' 



3(0fc - 4>k' 



sin^(f (fc-fc')) 



(25) 



The quasiperiodic boundary condition (j)N+i = f/*! + ^nw 
is employed. In fact, the simulation treats the phases as 
compact variables and the winding number w is com- 
puted by summing up the phase differences (between 
neighboring spins) along the chain. These phase differ- 
ences are defined modulo 2tt and we use the smaller angle 
- a procedure which should yield accurate results for suf- 
ficiently small time steps (see also comments in section 
HVBl) . 

The effective charging energy at gate charge zero can 
be computed from the winding numbers w of the paths 
(j) and the corresponding weights - the probabilities 
for a given w 



Ec 



27r2 E^w^P^ 



"G=0 



PEc Y.^P^ 



(26) 



Results for non-zero gate charge can be obtained from 
the expectation values (•••)«, in the different winding 
sectors w: 



[O] 



(27) 



Because of cancellation effects due to the phase factors 
in Eq. (|27p . the error bars can grow rapidly as a function 
of nc, especially at small a, where the P^-distribution is 
very broad. This manifestation of a sign problem restricts 
the parameter range for which accurate results can be 
obtained. 
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FIG. 2: Comparison of ^^^^ M{k, 0) defined in Eq. JgO]) 
(lines) and Monte Carlo results (symbols) for Xm (Eq. ((Sjl , 
= 0). The numerical data, which include error bars, were 
obtained for AtEc = 0.05. 
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FIG. 3: M{K,2nnG) defined in Eqs. ([20l) (lines) and 

Monte Carlo results (symbols) for Xm (0) and several val- 
ues of Q. The numerical data were obtained for l3Ec ~ 5 • 10'' 
and AtEc = 0.05; error bars at no = 0.125, 0.25, 0.375 and 
0.5 indicate the accuracy. 



We use WolfF cluster updates with efficient treatment 
of long range interactions [23lf23 | in our Monte Carlo 
simulations. This algorithm builds clusters in a time 
O(A^logA^) and allows us to simulate chains of up to 
10^ spins. For a discretization step AtEc = 0.05 
it is therefore possible to reach temperatures down to 
l3Ec = 5 • 10^. 



IV. RESULTS 

A. Magnetic susceptibility 

In Fig. [2] we compare the Monte Carlo results for the 
magnetic susceptibility for uq — against the scaling 
form dH), The plot for na = ^ ^ t^) looks pretty 
much the same as Fig. [21 Note that we applied relations 
([7]), ([5]), pH)) to find the value of the scahng parameter k 
corresponding to a given a and l3Ec ■ The agreement be- 
tween the analytical prediction and Monte Carlo results 
is within - 0.5% for PEc = 5 x 10", n = 1, 2, 3 and 
within - 1 - 2% for l3Ec = 5 x 10", n = 4, 5. Among 
other things, this indicates that the asymptotic series ([7]) 
truncated at fc = 1 approximates the Kondo temperature 
E* with an estimated error < 1% for any a > 1. This 
level of accuracy seems reasonable given that the k = 1 
term leads to a correction of approximately 10% for these 
a. 

As has been mentioned above there is no strong the- 
oretical argument that relation ((2^ is satisfied for 6 £ 
[0,7r]. However, as shown in the appendix, the saddle 
point approximation for is consistent with the hy- 
pothesis that the relation holds true for any 9 G [0,7r]. 
Further evidence that Eq. ([2^ is satisfied for all 9 comes 
from the study of the zero-temperature limit of the scal- 
ing function Fm{K,9). Namely, as follows from ([3]), (O, 



there exists a limit 

lim {KFra{K,e)) =/„(0) . (28) 

K — *OC 

Here fm (9) is some function of the single variable 6 which 
possesses a singular behavior as 6* — > tt. Using results of 
Refs. [21112^ one can show that 

fmiO)^^^ as Tr-9^+0. (29) 

At the same time, with the WKB approximation it is not 
difficult to prove the relation 

lim {kM [k, 9)) = -^ . (30) 
K^oo ^ ' sm(t^) 

Apparently Eqs. ([25)1 and ([50)1 are consistent with 
the conjecture Fm{n,9) — M{k,9). 

For the above reasons, we found it natural to test 
Eq. (|2^ numerically. The result is depicted in Fig. [31 
The missing data points for a = 1.2 and 1.3 fell victim 
to the sign problem mentioned in the previous section. 
This sign problem prevents us from investigating the sin- 
gularity at 9 = TT, which appears in the limit T ^ 0. 

B. Effective charging energy for zero gate charge 

First of all it is instructive to determine a range of 
temperatures appropriate for numerical tests of the scal- 
ing relation (jlip . The leading temperature-dependent 
correction to Eq. pTl) readily follows from the saddle 
point calculations presented in the appendix. Those 
calculations suggest (see formula (|A16P ) to replace the 
temperature-independent constant L = L{a) (|13p by 

L{a,[iEc)^L{a)+Li{^) , (31) 
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with 



Li{x) = 2e^ Ei(-a::) + 2 log{xe 



(32) 



For /3Ec — 5-10^, the term Li turns out to be fairly large 
(« 1.0). For this reason, it makes no sense to use Monte 
Carlo data with PEc — 5 • 10^ to test the relation pT|) . 
In the case (3Ec = 5 • 10^, Li « 0.2, so the correction 
is small but still visible and it is useful to take it into 
account. For PEc = 5 • 10" with n > 3, the term Li 
in Eq. (pij) is smaller than the error bars of our Monte 
Carlo data (< 0.03), and thus negligible. 

In Ref. [11] it was noted that the disagreement between 
previous Monte Carlo results [iill^ for E^ at a > I 
were due to lattice effects. Therefore it is very important 
to take proper care of the systematic errors introduced 
by the imaginary time discretization in the calculation of 

Eh- 

As a matter of fact, the strong effect of the discretiza- 
tion on the topological susceptibility is well known in 
the context of the 2D 0(3) nonhncar a-model (17ll8ill9i |. 
The primary reason of this effect is that strictly speaking 
the concept of winding number in a lattice formulation 
breaks down. In a discretized theory every field configu- 
ration can be continuously transformed into any other. If 
lattice configurations are sufficiently smooth, an unam- 
biguous topological charge may be assigned. Conversely, 
for field configurations containing large fluctuations the 
interpolation is not unique and the winding number defi- 
nition becomes ambiguous. For finite Ec the discretized 
Coulomb term in the lattice action ([25]) suppresses the 
lattice configurations containing large fluctuations. How- 
ever, it is possible to show that if a AtEq > 1 then start- 
ing with non-zero charged conflgurations one can contin- 
uously lower the lattice action Axv to zero by a local 
minimization in the spin variables [25|. For this reason 
the divergent constant L in Eq. pT|) (which is not uni- 
versal and depends on details of the discretization of the 
action (21)) is a singular function of the lattice parameter 



S = a AtEc 



at (5 = 1. Namely, 



i cx -2 log(l - S) 



(33) 



(34) 



In Fig.m 2Tr^E' "^ith n-G = is plotted as a function 
of AtEc for l3Ec — 5 • 10"^ and some values of a. An 
analysis of the plotted Monte Carla data suggests that 
the effect of the lattice discretization for (5 <C 1 can be 
accounted for by means of the following modiflcation of 
Eq.®: 

L{a,PEc,S) = L{a) + Li{j^) 

-2 log(l - (5) + 2(7(5 . (35) 



Here L{a) and Li{x) are given by Eqs. (1131) ^'^d (l32l) . re- 
spectively, and C is a fitting parameter. The continuous 
lines in Fig. 21 correspond to C = 2.8. 
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FIG. 4: Illustration of lattice effects for 
Monte Carlo data correspond to a 
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{no = 0). The 
1.6 and PEc = 



5 • 10 . The lines show the analytical results corresponding to 
the lattice fitting parameter C = 2.8 in Eq. (I35|l . 



In order to minimize the impact of the non-universal 
temperature and lattice effects discussed above, we first 
study the quantity (jig — 0) 



fnm{a) 



E*r 



2Tr^E* 



E, 



c 



/3Ec=5-10" 



2Tr^E* 



(36) 



l3Ec=5-10" 



As follows from Eq. (fTl]) . the anomalous terms cancel out 
for such differences. Hence one should expect that the 
functions /„,„ possess the normal scaling behavior of the 
form 



/nm(a) Ft{Kn,0) - i^t(K,„,0) ( n, TO, 



a 



oo 



(37) 



with Kn^m the values of the scaling parameter n pop cor- 
responding to a given a and PEq = 5 • lO"'™. Also 
notice that the non-universal temperature corrections to 
the scaling behavior (|37l) are expected to be sufficiently 
small provided n, m > 3. Since the most accurate Monte 
Carlo data are those for n = 3, we plot in Fig.Elthe result 
for /43, obtained for a discretization step AtEc = 0.05, 
and the analytical prediction (Eqs. (|37l) . (HH)). The agree- 
ment is as good as it can be expected, given the errors 
on the Monte Carlo data, and a similar result is found 

for /32- 

In Fig.inithe prediction (HIl), (UHl) for = and L 
understood as in Eq. (|551) . is compared to the Monte 
Carlo data for different values of the parameters /SEc 
and AtEc- The value C = 2.8 has been used for all an- 
alytical curves. The agreement between the theoretical 
predictions and the Monte Carlo results shows that the 
scaling formulas (fTTl) . (fT3l) . (fT8l) for the effective charg- 
ing energy are correct. In particular, this implies that 
the leading power in the pre-exponential factor /(a) (see 
Eqs. (ESI) and ^) is . 
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FIG. 5: Comparison of the Monte Carlo data and analytical 
prediction for /43 defined in Eq. (I36|l . The numerical data 
correspond to AtEc ~ 0.05 and the analytical result was 
obtained using Eqs. (P7|l . HHJ. 



Carlo data and the theoretical predictions over a whole 
range of temperatures (5 • 10^ < /3i?c < 5 • 10^), lat- 
tice spacings (0.005 < AtEc < 0.05) and dissipation 
strengths (1.2 < a < 2.4). These results essentially ver- 
ify the scaling formulas for the effective charging energy 
for no = presented in Ref. [l^l, and thus also settle 
the controversy surrounding the pre-exponential factor 
/(a) in Eq. pS)) . Unfortunately the quality of the Monte 
Carlo data for E^{nG) with no 7^ was not sufficient 
for a meaningful test of the 6'-dependence of the scaling 
function Ft{K,6) pT|) against the theoretical prediction 

m- 

As for the magnetic susceptibility, our investigation 
suggests that some of the results from Refs. fi\w\ for 
no = and no = 1/2 can be generalized to all values of 
the gate charge. 
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APPENDIX A: SADDLE POINT 
APPROXIMATION 

In this appendix we consider the saddle point approx- 
imation for the model ^ in the limit Ec — > cxd. In the 
absence of the first term, the local minima of the Eu- 
clidean action A[(j)] in the sectors with winding number 
w > are saturated by the instanton solutions 20ll27t 



with z = e2'^^/'^ 



(Al) 



V. SUMMARY 

We have presented scaling formulas for the single elec- 
tron box in the limit of almost perfect transmission, as 
well as a conjecture for the gate charge dependence of 
the magnetic susceptibility. These analytical predictions 
have been tested by means of Monte Carlo simulations 
of the discretized action (corresponding to an XY spin 
chain with long range couplings) with up to 10^ spins. 
Because of the lattice effects, it was necessary to intro- 
duce a constant, C — 2.8, in order to compare the analyt- 
ical predictions to the simulation results for the effective 
charging energy. But this single fitting parameter was 
enough to obtain a good agreement between the Monte 



The set of complex parameters {z^ ■ l^ml < l}m=i ^^^d 
the real constant < 0o < 27r are instanton moduli. Let 
us introduce an external magnetic field in the system: 

1-13 

Ah[(f>] = A[(f>] + hV^ / dr cos (0(r)) . (A2) 
Jo 

We assume that h is sufhciently small so that the addi- 
tional term has no appreciable effect on the instanton so- 
lutions. Calculating the action .A;i[0] for the paths (|A1|) 
yields 

Ah[(t>inst] = iOw + an"^ w + ^ I3hy/a {-ly X 



n 



(A3) 
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where 9 is the topological angle ([5]). Within the saddle 
point approximation the partition function is given by 



oo ^ 
tu=0"''" 



(A4) 



where stands for the integral over the domain < 
4>o < 2tt, \zm\ < 1 (m = 1, . . .w). The integration mea- 
sure d^Aw should be computed as usual by integrating 
out the Gaussian fluctuations around the instanton solu- 
tions. It was calculated in Rcfs. [6.,2&]: 



dA^, 



Here 



(/3A) 



(27r2a)5+'" X 



, „ (A5) 

det II G(z„z,) II . 



m— 1 



Gizz,zj) = — 



(A6) 



and A is the perturbative ultraviolet cut-off ([6]). 

Expression (jA4p can be evaluated in closed form. In- 
deed, let us expand the exponent e~'^''['^'"='l into a power 
series of ph and then integrate explicitly over (f>Q. This 
yields 



■^inst(^) 



p=0 

with K — (3 Aae""'^ and 



z,w = i+E^ n / 



(A8) 



n {z^z^mf det||G(z„z;)|| 



27ri 



It is easy to see that Zp{X} can be interpreted as a (nor- 
malized) fermionic partition function: 



J VcDc Q-(ZG-^ c)-\(c,c) 

J VcDc e~(eG-ic) 
det(G-i + A) 



(A9) 



det(G-i) 



det(l + AG) 



where the operator G acting on the Grassman fields is 
defined by the equation 



(Gc)(z,z*) 



(AlO) 



dC A dC* 



ICKi 



27ri 



[ccY G(z,c)c(c,r) 



Hence, 



logZp = Splog(l + AG) = - ^ Sp(G^), (All) 



with 



Sp(G^) = 
[(1 



n 



A dz 



>l<i 



27ri 



{Zm ZmY 



Zizl){l - Z2Z3) ... (1 - ZsZl) _ 



(A12) 



The small instanton problem is explicit now. Indeed, 
equation (|A12p for s = 1 can be written in the form 



Sp(G) 



ki|< 



dzi A dz* 
I 27ri(l — zizl) 



-7£-V'(l+p),(Al3) 



where tp{u) = ^ logr(u) and the integral diverges as 
1^1 1 1. This limit corresponds to small instantons 
- when |zi| is close to one, the solution 0inst (|A1[) 
with w = 1 is almost constant everywhere in the seg- 
ment [0, /3] except for a small neighborhood of the point 
T = ^arg(zi) where the instanton is localized. Any 
regularization (denoted by j' in (|A13p ) is essentially an 
instruction to suppress contributions from the small-size 
instantons. As a matter of fact, the Coulomb term in 
the action ([2|) provides such a suppression for any finite 
f3Ec » 1. Thus the integral in (|A13p should be under- 
stood as M - 



dzi A dzl e 



|zi|<l 



27ri 



1 



ZlZ, 



Here 



A 



(1) _ 



1 + zizl 
PEc l-zizl 



(A14) 



(A15) 



is the value of the Coulomb term calculated on the one 
instanton solution (w = 1). Then 

Sp(G) = 1 [L + Li + login) - + p), (A16) 

where L = 27r^a-|-0(loga) is some divergent temperature 
-independent constant. The function Li(x) is given by 
Eq. p2|) . It vanishes in the limit /3Ec' — > 00. 

All integrals (jA12[) with s > 1 are finite. Integrating 
over the phases of Zm in (IA12P first, one obtains [xm = 

\z„r), 



Sp(G«) 



n 



E 

n=0 



dXn 



dx 



1 — Xl 



Cp(s) , 



(A17) 



where C.p{s) — '^'^^oi^+p)^'^ is the generalized Riemann 
^-function. Finally, .^inst(^) can be expressed in terms of 
the conventional Bessel function, 

^i„.t(/i) = e5i^ (_|^^^/^(^/,^),(A18) 
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with A = Ace"'^. 

Some explanations are in order at this point. First, the 
corrections to the saddle point formula (|A4p are small 
in a^^ but they contain ultraviolet divergent integrals. 
Taking these divergences properly into account leads to 
a renormalization of the bare coupling go = (27r^a)^^ 
in Eq. (jA18|) . The renormalization trades go for the 
"running coupling constant" g{K) subject to the RG flow 
equation 



d5 



2g'+Ag-' + 0ig^) 



(A19) 



Second, in the above derivation we took into account 
the instantons with positive winding numbers only. One 
can indeed ignore the "anti-instanton" configurations 
corresponding to w < provided 9m(0) ^ 1. Obviously, 
the contribution of anti-instanton fluctuations dominates 
when Qm(— 0) 1. In this case the saddle point approx- 
imation leads to (jAlSp with A = kc'^. Both limiting cases 
are neatly captured by the same formula if A is under- 
stood as 



X — 2k cos 9 



(A20) 



Equation (|A20p is not easy to justify on general grounds 
for 9^1, but can be supported by means of the exact 
result in the case h = 0. 

According to Ref. [l^l , the scaling behavior of the par- 
tition function ([T]) looks as follows 



^scal — e 



Ln cos 6 



Z{k,9) {/3Ec, 



,(A21) 



where L and Z are given by Eqs. (|T3|) and (|T9|) . respec- 
tively. We consider now the high-temperature (small-^t) 
expansion of (|A21[) . Here it is useful to introduce the 
running coupling constant g = g{K) as a solution of the 
equation 



-1 - — 

e 2g 



(A22) 



Note that g ^ 1 and it solves the RG flow equation (|A19[) 
within the two-loop approximation. After a change of 
variables 

y^gx-logg, (A23) 
the differential equation (jlSp can be written in the form 

-dl + c^ + {gKC0s9f + SU{x)] ^-(0;) = 0, (A24) 



where 



SU (x) — exp 



(A25) 



For ~ 1 the term 6U is 0{g), so the solution "ii^ is 
approximated by the Macdonald function 

^+ix) = 2g^ K2gncosei2e^) (l + 0(g)) . (A26) 
The normalization factor in front of K is fixed by match- 
ing (jA26[) to the asymptotic condition (|17p . As (— .t) ^ 
1, takes the form 



The coefficient C_ is simply related to the Wronskian in 
Eq. HI]), namely 

= 2k cos 6*3-"'=°"^ C_ . (A28) 



In view of (|A24p . one can show that the high-temperature 
expansion of Zgcai (|A21[) has the form: 



•^scal — ^ 



h:(L— log g) cos 9 



V9 



(A29) 



m=0 



where the coefficients Z„i{k,9) admit systematic expan- 
sions in powers of k. In particular 



Zo{k,9) = 



r(l + 2k cos ( 



(A30) 



It is straightforward to check that the saddle point par- 
tition function Zinst(O) is in agreement with (|A30[) pro- 
vided A in equation (jAlSp is given by (|A20p and the effect 
of the renormalization of go = (27r^Q!)^^ is taken into ac- 
count. 

Let us return to Eqs. (jAlSp . (|A20p to find the magnetic 
susceptibility within the saddle point approximation: 



(A31) 



277% g{K) (1 + 2k cos 6*) ' 

Again the effect of renormalization has been accounted 
for in (jA3ip . It is easy to see that formula (|A3ip is 
in agreement with the hypothesis that relation is 
satisfied for any 9. Indeed, as follows from the above 
analysis, in the regime k ^ 1 the potential U{y) — 
K^ exp ( e^ ) — K^ sin^ (9) in the differential equation ([15]) 
has the effect of a rigid wall at some point yo ~ ~ log5i 
i.e., to the left from this point, the potential is negligi- 
ble, but to the right, it grows very fast, so the solution 
^'+(t/) is essentially zero. More precisely, when y is be- 
low but close to yo the solution is approximated by 
a linear function, '^+{y) ~ g^^ {yo — y). Within this 
"rigid wall" approximation A+Ijjq) « 0, 9j,A+(yo) ~ 
and s(yo) > 1- Therefore, the RHS of Eq. ((201), calcu- 
lated at 2/ = yo, can be approximated by A-{yo). In the 
calculation of the solution needed for the in- 

tegral (pij) can be replaced by its asymptotics ([TB)l . This 
yields the relation 



M(k, 9) W 277^0! X 



(inst) 



(A32) 



*+(a;) ^ C+ e^'=9™'^» -I- C_ 6"=^'^^''=°''^ . (A27) where x^"'*^ is given by (IMD) 
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